FTUAM/08-18 
IFT-UAM/CSIC-08-58 



Lepton flavour violation in constrained MSSM-seesaw models 



E. Arganda. M.J. Herrero 

Departamento de Fisica Teorica C-XI and Instituto de Fisica Teorica C-XVI, Universidad Autonoma 

de Madrid, Cantoblanco, E-28049 Madrid, Spain 

00 '. 

^— «j \ We calculate the predictions for lepton flavour violating (LFV) tan and muon decays, Ij — > h^, 

^ . Ij 3li, n ~ e conversion in nuclei and LFV semileptonic tan decays r — * jj,PP with 

^ ■ PP = Tv+-K-,n°n°,K+K-,K°K° t ^ fiP with P = 7r°,7?,V and r ^ /xV with V = p°,(j), 

performing the hadronisation of quark bilinears within the chiral framework. We work within 
O ' a SUSY-seesaw context where the particle content of the Minimal Supersymmetric Standard 

' Model is extended by three right-handed neutrinos plus their corresponding SUSY partners, 

and where a seesaw mechanism for neutrino mass generation is implemented. Two different 
scenarios with either universal or non-universal soft supersymmetry breaking Higgs masses 
at the gauge coupling unification scale are considered. After comparing the predictions with 
■ present experimental bounds and future sensitivities, the most promising processes are par- 

l' ticularly emphasised. 

Oh' 

1 LFV within SUSY-seesaw models 

. The current knowlegde of neutrino mass differences and mixing angles clearly indicates that 

\^ \ lepton flavour number is not a conserved quantum number in Nature. However, the lepton 

flavour violation (LFV) has so far been observed only in the neutrino sector. One challenging 
task for the present and future experiments will then be to test if there is or there is not LFV 
in the charged lepton sector as well. 



o 
d 



OO , Here we focus in the Minimal Supersymmetric Standard Model (MSSM) enlarged by three 

' right-handed neutrinos and their SUSY partners where potentially observable LFV effects in the 

charged lepton sector are expected to occur. We further assume a seesaw mechanism for neutrino 
^ ■ mass generation and use, in particular, the parameterisation proposed in^ where the solution 

c3 ' to the seesaw equation is written as m^i = YpV2= \Jrr^^ R\frr^^ulj[^^. Here, R is defined 

by Oi {i = 1,2,3); Ui(2) = cos(sin)/3, v = 174 GeV; m^'^^ = diag (mi,^ , mj^j , JTi^^g ) denotes the 

three light neutrino masses, and m^^^ = diag (rriAr-^ , m7V2 , TVs ) the three heavy ones. ?7mns is 
given by the three (light) neutrino mixing angles ^12,^23 and ^13, and three phases, 5,4)1 and 
<^2- With this parameterisation is easy to accommodate the neutrino data, while leaving room 
for extra neutrino mixings (from the right-handed sector). It further allows for large Yukawa 
couplings ~ 0(1) by choosing large entries in m'^^^ and/or di. 

The predictions in the following are for two different constrained MSSM-seesaw scenarios, 
with universal and non-universal Higgs soft masses and with respective parameters (in addi- 
tion to the previous neutrino sector parameters): 1) CMSSM-seesaw: Mq, Mi/2j ^0 tan/?, 
and sign(;u), and 2) NUHM-seesaw: Mq, ^0 tan/3, sign(//), Mh^ = Mo(l + 5i)^/^ and 

= Mq{1 + (^2)^/^. All the predictions presented here include the full set of SUSY one-loop 
contributing diagrams and we do not use the Leading Logarithmic (LLog) nor the mass insertion 
approximations. The hadronisation of quark bilinears is performed within the chiral framework, 
using xPT and R^T. This is a very short summary of several publications ^ ^ ^ to which we 
refer the reader for more details. 
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Figure 1: t ^ fi^ and CR(^ — e, Ti) as a function of mjvj. The predictions for SPS la (dots), lb (crosses), 
2 (asterisks), 3 (triangles), 4 (circles) and 5 (times) are included. On the upper horizontal axis we display the 
associated value of {Y^)^^. In each case, we set On — 5°, and 9i — 0. The upper (lower) horizontal line denotes 

the present experimental bound (future sensitivity). 



2 Results and Discussion 

We focus on the dependence on the most relevant parameters which, for the case of hierarchical 
(degenerate) heavy neutrinos, are: the neutrino mass ("t-at), tan/3, 9i and 62- We also 
study the sensitivity of the BRs to ^13. The other input seesaw parameters m-ATj, and 
^3, play a secondary role since the BRs do not strongly depend on them. The light neutrino 
parameters are fixed to: m^^ = ^"^soi + ''^ui^ "^-^3 = ^"^atm + "^^ii ^"^soi = 8 x 10~^eV^, 
Aml^^ = 2.5 X 10-3 eV^ m^^ = W'^ eV, 6*12 = 30°, 623 = 45°, 6*13 < 10° and 6 = (pi = (p2 = 0. 

The results for the CMSSM-seesaw scenario are collected in Figs. 1 through 5. In Fig. 1, 
we display the predictions of BR(r — fij) and CR(^ — e, Ti) as a function of the heaviest 
neutrino mass m]\f,^ for the various SPS points, and for the particular choice 6i = {i = 1,2,3) 
and ^13 = 5°. We have also considered the case of degenerate heavy neutrino spectra (not 
shown here). In both scenarios for degenerate and hierarchical heavy neutrinos, we find a strong 
dependence on the the heavy neutrino masses, with the expected behaviour |m7v logm^vp of the 
LLog approximation, except for SPS 5 point, which fails by a factor of ~ 10^. The rates for the 
various SPS points exhibit the following hierarchy, BR4 > BRib > BRia > BR3 > BR2 > BR5. 
This behaviour can be understood in terms of the growth of the BRs with tan (3, and from the 
different mass spectra associated with each point. Most of the studied processes reach their 
experimental limit at mj\^^ € [10^^,10^^] which corresponds to Y^^'^'^ ~ 0.1 — 1. At present, 
the most restrictive one is /x — 67 (which sets bounds for SPS la oi ruN^ < W^^ — 10^^ GeV), 
although fi — e conversion will be the best one in future, with a sensitivity to mTVg > 10^^ GeV. 

Fig. 2 shows the behaviour of the six considered LFV r and decays, for SPS 4 point, as 
a function of \0i\, for various values of arg^i. We see clearly that the BRs for < |^i| < vr and 
< arg^i < 7r/2 can increase up to a factor 10^ — 10^ with respect to 6i = 0. Similar results have 
been found for 62, while BRs are nearly constant with 9^ in the case of hierarchical neutrinos. 
The behaviour of CR(/i — e, Ti) with 6i is very similar to that of BR(// — s- ej) and BR(/i 3e). 
For instance. Fig. 3 shows the dependence of CR(/i — e, Ti) with 62, and illustrates that for large 
62, rates up to a factor ~ 10^ larger than in the 9i = case can be obtained. 

In Fig. 4 we show the dependence of /i — > 67, ^ 3e and /i — e conversion on the light 
neutrino mixing angle ^13. These figures clearly manifest the very strong sensitivity of their 
rates to the ^13 mixing angle for hierarchical heavy neutrinos. Indeed, varying ^13 from to 10° 
leads to an increase in the rates by as much as five orders of magnitude. 

On the other hand, since ;U ^ 67 is very sensitive to ^13, but BR(t ^7) is clearly not, 
and since both BRs display the same approximate behaviour with rriN.^ and tan (5, one can study 
the impact that a potential future measurement of ^13 and these two rates can have on the 




Figure 2: Dependence of LFV r and ^ decays with l^i] for SPS 4 case with arg(Si) = 0, vr/lO, 7r/8, 7r/6, 7r/4 in 
radians (lower to upper lines), (m]Vi , mjv^ , mjv3) = (10*, 2 x 10*, 10^*) GeV, 62 — 63 = 0, 613 = and m,^-^ = 0. 

The horizontal lines are the present experimental bounds. 
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Figure 3: CK{fi — e, Ti) as a function of I6I2I, for SPS la case with arg6l2 ~ {0, 7r/8, 7r/4, Stt/S, 7r/2} (dots, 
crosses, asterisks, triangles and circles, respectively), mjVi = (10^^°, lO'^'^, lO'^^) GeV, 613 = 5°. The upper (lower) 
horizontal line denotes the present experimental bound (future sensitivity). 



knowledge of the otherwise unreacheable heavy neutrino parameters. The correlation of these 
two observables as a function ofm]\[^, is shown in Fig. 5 for SPS la. Comparing these predictions 
for the shaded areas along the expected diagonal "corridor", with the allowed experimental 
region, allows to conclude about the impact of a ^13 measurement on the allowed/excluded m]\f^ 
values. The most important conclusion from Fig. 5 is that for SPS la, and for the parameter 
space defined in the caption, an hypothetical ^13 measurement larger than 1°, together with 
the present experimental bound on the BR(^ — > 67), will have the impact of excluding values 
of rriTVa ^ 10^"^ GeV. Moreover, with the planned MEG sensitivity, the same ^13 measurement 
could further exclude rriN^ > 3 x 10^^ GeV. 

The numerical results for the NUHM-seesaw scenario as a function of Mq = M1/2 = MsusY 
are collected in Figs. 6 and 7. The behaviour of the predicted rrijjo as a function of Mgugy is 
shown in Fig. 6 (left panel). The most interesting solutions with important phenomenological 
implications are found for negative 5i and positive 82- Notice that, for all the explored 81^2 
values, we find a value of mjjo that is significantly smaller than in the universal case {61^2 = 0). 

In Fig. 6 (right panel) the various contributions from the 7-, Z-, Higgs mediated penguins 
and box diagrams as a function of MsusY are shown. Here, we choose 5i = —1.8 and 82 = 0. We 
observe a very distinct behaviour with MsusY of the Higgs-mediated contributions compared 
to those of the CMSSM case. In fact, the Higgs-mediated contribution can equal, or even 
exceed that of the photon, dominating the total conversion rate in the large Mq = Mi ^2 region. 



t 10" 

3 



rrifj = (10^°,10",10^* ) GeV 
m^i = 10"^ eV 
01 = 



SPS 1a • 
SPS lb 

SPS 2 * 

SPS 3 

SPS 4 

SPS 5 



n 



10-3 








rrifj = {10"',10",10 






io'° 


m^,i = 10"^ eV 
= ^ 




000000000000000=000 . 


10-1^ 
10"''^ 














10"''* 






SPS la • 

SPS lb + 


10-^^ 






SPS 2 >K 
SPS 3 » 


10-1^ 






SPS4 o 
SPS 5 


10-'^ 
io'^ 








2 


4 

9l3 


6 8 10 

(°) 



m^, - (10 , 10 , 10 ) GeV 
6| - 

,0000000001 



.oQ°' 



lOOOO' 



lOOOOf 



oooooooooooooooooo 



!f»»" SPS la • 

SPS lb + 
SPS 2 « 
SPS 3 
SPS 4 o 
SPSS 



8l3 (°) 



Figure 4: BR(/i 67), BR(/i — > 3e) and CR(/i — e, Ti) as a function of 813 (in degrees), for SPS la (dots), lb 
(crosses), 2 (asterisks), 3 (triangles), 4 (circles) and 5 (times), with 6i = and mjv. = (10^°, 10^^, 10^*) GeV. The 
upper (lower) horizontal line denotes the present experimental bound (future sensitivity). 




Figure 5: Correlation between BR(/i -^67) and BR(r fi'y) as a function of miv.,, for SPS la, and impact 
of 9is. The areas displayed represent the scan over 9i. From bottom to top, the coloured regions correspond to 
6'i3 — 1°, 3°, 5° and 10° (red, green, blue and pink, respectively). Horizontal and vertical dashed (dotted) lines 

denote the experimental bounds (future sensitivities). 



These larger Higgs contributions are the consequence of their exclusive SUSY non-decoupling 
behaviour for large MsusY; and of the lighter Higgs boson mass values encountered in this 
region, as previously illustrated in Fig. 6. 

In Fig. 7 we display the predicted /u — e conversion rates for other nuclei, concretely Al, Ti, 
Sr, Sb, Au and Pb, as a function of MsusY- We clearly see that CR(/U — e, Sb) > CR(/u — 
e, Sr) > CR(^ - e, Ti) > CR{n - e, Au) > CR(// - e, Pb) > CR{n - e, Al). The most 
important conclusion from Fig. 7 is that we have found predictions for Gold nuclei which, for 
the input parameters in this plot, are above its present experimental bound throughout the 
explored MsusY interval. Finally, althought not shown here for shortness, we have also found 
an interesting loss of correlation between the predicted CR(/i — e, Ti) and BR(// — >■ 67) in the 
NUHM-seesaw scenario compared to the universal case where these are known to be strongly 
correlated. This loss of correlation occurs when the Higgs-contributions dominate the photon- 
contributions and could be tested if the announced future sensitivities in these quantities are 
reached. 

The corresponding predictions for 02 = 2.9e*'^/^ of the nine LFV semileptonic r decays 
studied in this work as a function of MgusY are shown in Fig. 8. In this case, we work with 
61 = —2.4 and 82 = 0.2, that drive us to Higgs boson masses around 150 GeV even for heavy 
SUSY spectra. In this Fig. 8 we can see that, the choice of 62 increase all the rates about two 
orders of magnitude respect to the case 6i = 0, not shown here for brevitiy. BR(t /ivr^vr") 
and BR(t — > up) get the largest rates and, indeed, the predictions of these two latter channels 
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Figure 6: Left panel: Mass of mjjo as a function of AIq ~ M1/2, for fixed values of Si — { — 1.8, —1.6, —1, 0} 
(respectively crosses, asterisks, triangles and circles), with mjv. = (10^", 10^^, 10^*) GeV, 6i = 0, Ao = 0, tan/3 = 

50 and Ora = 5°. Right panel: Contributions to CK{fi — e, Ti): total (dots), 7-penguins (diamonds), .Z-penguins 
(asterisks), Jf-penguins (crosses) and box diagrams (times) as a function of Mo(= M1/2) for the NUHM case with 

51 = -1.8, S2 = 0, tan/3 = 50, mjv, = (10^°, lO", 10^*) GeV, 613 = 5° and .R = 1 {Oi = 0). The upper (lower) 

horizontal line denotes the present experimental bound (future sensitivity). 
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Figure 7: fj. ~ e conversion rates as a function of Afo = M1/2 in the NUHM-seesaw for various nuclei: Sb, Sr, 
Ti, Au, Pb and Al nuclei (diamonds, triangles, dots, asterisks, times and crosses, respectively), with mjv. = 



(10'°,10",10''*) GeV, Ao = 0, tan/3 = 50, 613 = 5° 



0, Si — —1.8 and S2 — 0. From top to bottom, the 



horizontal dashed lines denote the present experimental bounds for CR(/^ — e, Ti) and CR(/i — e, Au). 



reach their present experimental sensitivities at the low -MsusY region, below 200 GeV and 250 
GeV respectively, for this particular choice of input parameters. 

In Fig. 9 we plot finally the predictions for BR(r — > ijlK~^K~) and BR(r — > ixrj) as a 
function of one the most relevant parameters for these Higgs-mediated processes which is the 
corresponding Higgs boson mass. 

Firstly, we see that the approximate (see the approximate formulae in^S) and exact results of 
the Higgs contribution agree within a factor of two for both channels, but the agreement of the 
full result with respect to the Higgs contribution is clearly worse in the case of r — > fiK~^K~ than 
in r — > iir]. In the latter, the agreement is quite good because the Z-mediated contribution is 
negligible, and this holds for all MsusY values in the studied interval, 250 GeV < -MsusY < 750 
GeV . In the first, it is only for large MsusY that the .ff-mediated contribution competes with 
the 7-mediated one and the Higgs rates approach the total rates. For instance, the predictions 
for BR(r — > fiK~^K~) shows that for -MsusY = 750 GeV and m^o = 160 GeV the total rate is 
about a factor 2 above the Higgs rate, but for rrifjo = 240 GeV it is already more than a factor 
5 above. 

In this figure we have also explored larger values of mjvg and tan /?, by using in those cases the 
approximate formula, and in order to conclude about the values that predict rates comparable 
with the present experimental sensitivity. We can conclude then that, at present, it is certainly 
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Figure 8: Predictions of BR(r fJ'PP) and BR(r /iP) as a function of Msusy in the NUHM scenario for a 

large r — /i mixing driven by 62 = 2.9e"^/''. 
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Figure 9: Predictions for BR(r fj,K^ K ) and BR(t fir]) as a function of rrijjo in the NUHM scenario. 



T ^ fir] the most competitive LFV semileptonic tau decay channel. The paremeter values that 
provide rates being comparable to the present sensitivities in this channel are tan /? = 60 and 
m^Vs = 10^^ GeV which correspond to {6^21 ~ 2. 

Interestingly, the most competitive channels to explore simultaneously LFV r — ^ transitions 
and the Higgs sector are r — > fir], r — >• iirj' and also r — fiK^K" . Otherwise, the golden channels 
to tackle the Higgs sector are undoubtly t ^ firj and r fir]'. On the other hand, the rest of 
the studied semileptonic channels, r /i7r"*"7r~, etc., will not provide additional information on 
LFV with respect to that provided by r ^ ^7. 

In conclusion, we believe that a joint measurement of the LFV branching ratios, the ji — e 
conversion rates, ^13 and the SUSY spectrum will be a powerful tool for shedding some light 
on the otherwise unreachable heavy neutrino parameters. Futhermore, in the case of a NUHM 
scenario, it may also provide interesting information on the Higgs sector. It is clear from this 
study that the connection between LFV and neutrino physics will play a relevant role for the 
searches of new physics beyond the SM. 
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